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Stimulated by recent experimental discoveries, triaxial strongly deformed (TSD) states in 158 Er 
at ultrahigh spins have been studied by means of the Skyrme-Hartree-Fock model and the tilted- 
axis-cranking method. Restricting the rotational axis to one of the principal axes - as done in 
previous cranking calculations - two well-defined TSD minima in the total Routhian surface are 
found for a given configuration: one with positive and another with negative triaxial deformation 
7. By allowing the rotational axis to change direction, the higher-energy minimum is shown to be 
a saddle point. This resolves the long-standing question of the physical interpretation of the two 
triaxial minima at a very similar quadrupole shape obtained in the principal axis cranking approach. 
Several TSD configurations have been predicted, including a highly deformed band expected to cross 
lesser elongated TSD bands at the highest spins. Rs transitional quadrupole moment Q t ~ 10.5 eb 
is close to the measured value of ~lleb; hence, it is a candidate for the structure observed in 
experiment. 



PACS numbers: 21.60.Jz, 21.10.Re, 21.10.Ky, 27.70.+q 



While the majority of nuclei have axially symmetric 
shapes, evidence for triaxial nuclear deformations has 
been elusive. The clearest signatures come from the 
gamma-ray spectroscopy of rotating nuclei. The defor- 
mation of a quantum object, such as molecule or atomic 
nucleus, enables the system to specify an orientation. 
The quantized motion of this degree of freedom gener- 
ates the sequences of rotational levels - the rotational 
bands [3, HI • If the system is triaxial, the associated ro- 
tational bands show specific features that allow for distin- 
guishing it from an axial one. In the case of nuclei, the ap- 
pearance of the wobbling [H, H, 01 and spin-chirality 0, [1| 
rotational modes are experimental signatures of triaxial- 
ity 

Triaxial shapes are expected to appear more frequently 
at high spin because of the tendency of aligned high-j 
quasiparticles to drive rotating nuclei towards triaxiality 
due to their spatial density distributions d, Q • In ad- 
dition, pairing correlations - which generally favor more 
symmetric shapes - are quenched at high spins and en- 
hance the high-j alignment effect Q . Consequently, with 
increasing spin, nuclei are predicted to go through non- 
axial shapes before they eventually fission (see, e.g., 0). 

Recent experiments fToj - TTsI ] have reached ultrahigh 
spins of about 65h in nuclei around 158 Er. R has been 
observed that with increasing angular momentum, the 
rotational bands terminate and nuclei assume weakly de- 
formed oblate shapes, as evidenced by the irregular level 
spacings. At ultrahigh spins they return to collective 



rotation characterized by regular rotational bands, con- 
sistent with the early prediction (T^). Cranked Nilsson- 
Strutinsky (CNS) calculations suggest that the observed 
bands in 158 Er are based on one of the three triaxial 
strongly deformed (TSD) minima in the potential-energy 
surface (see Ref. [16| and references cited therein). For 
the lowest minimum, TSD1, the calculated value of the 
transitional quadrupole moment Q t ~ 7.5 cb consider- 
ably underestimates the observed value of ~ 1 1 eb [r| . 
This has led to the suggestion that the observed band in 
158 Er may be associated with either the minimum TSD2, 
which has a similar quadrupole deformation parameter £2 
as TSD1 but negative 7, or with the band TSD3, which 
has a larger triaxial deformation p~6j | . 

Most of the existing high-spin calculations in the mass- 
160 region assume that the axis of rotation coincides with 
one of the principal axes of the triaxial potential, which is 
commonly referred to as principal-axis cranking (PAC). 
The CNS calculations use the microscopic-macroscopic 
method, which combines a shell correction derived from 
a phcnomenological potential with the deformation en- 
ergy of a rotating liquid drop [l7l - [l9j . It is common to 
choose the x-axis as the rotational axis and let the tri- 
axiality parameter cover the range —120° < 7 < 60°. 
In the Lund convention, which we adopt in this Letter, 
the three sectors [-120°, -60°], [-60°, 0°], and [0°,60°] 
represent the same triaxial shapes but represent rotation 
about the long, medium and short axis, respectively. The 
TSD1 and TSD2 minima in the CNS calculations [r| cor- 
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respond to similar values of £2 and | — ^ | which means that 
their shapes are nearly the same. The opposite sign of 7 
means that TSDl rotates about the short axis and TSD2 
about the medium axis. This raises the question of their 
physical interpretation, e.g., whether the higher of the 
two minima obtained in PAC is stable with respect to a 
reorientation of the rotational axis. 

In this Letter we address this question by means of the 
tilted axis cranking (TAC) method [I(|[2l[, which con- 
siders the general orientation of the axis of rotation with 
respect to the principal axes of the nuclear quadrupole 
moment. We investigate the structure of ultrahigh-spin 
TSD minima in 158 Er by using two approaches: the 
shell-correction tilted axis cranking method (SCTAC) 
[20| . which is based on the phenomenological Nilsson 
potential, and, for the first time, the three-dimensional 
self-consistent Skyrmc-Hartree-Fock version of tilted axis 
cranking (SHFTAC) developed in Ref. [22|. Employing 
a fully self-consistent rotating mean field - including the 
full rotational response due to cranking - is expected 
to improve the reliability of calculations in the realm of 
ultrahigh-spin states. 

SHFTAC is based on the symmetry-unrestricted solver 
HFODD (v2.49s) [23|, which has been successfully applied 
to the description of chiral bands in 132 La [23]. In the 
particle-hole channel, we use the Skyrme energy den- 
sity functionals SkM* Q and SLy4 [Hj], the latter of 
which has been supplemented with Landau parameters 
(SLy4L) 0, [27]] ■ The total energy of the system £ is 
obtained by integrating the total energy density over 
spatial coordinates. We have used 1,000 deformed har- 
monic oscillator basis states with haj± — 10.080 MeV and 
fiwy = 7.418 MeV. At ultrahigh spins, pairing is negligi- 
ble; hence, it has been ignored in SCTAC and SHFTAC. 
As discussed in earlier SHFPAC (a PAC limit of SHF- 
TAC) calculations 0, quadrupole polarization at high 
spin - both axial and triaxial - is very well described by 
unpaired theory The SCTAC calculations give a strong 





FIG. 1: (Color online) Schematic picture of a TSD shape. 
The angles 9 (between the :r-axis and the rotational axis) and 
a (between u> and J) are defined in the x — y plane. The 
short, medium, and long axes are denoted by x, y, and z, 
respectively; that is, the plotted shape corresponds to 7 > 0. 

increase of the Routhian £ w = £ — ui ■ J when the ro- 



tational axis is tilted toward the z-axis. For this reason 
we restrict the numerically extensive SHFTAC calcula- 
tion to the x — y plane spanned by the short and medium 
axes. The tilt angle 9 of the rotational axis is measured 
with respect to the short axis, as illustrated in Fig. [TJ 
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FIG. 2: (Color online) Lowest Routhians in 158 Er calculated 
by means of SCTAC (Nilsson potential) at fixed £2 and uj 
as a function of 7 for rotation about the x-axis (solid line) 
compared to those calculated at 7 = 15° and 22° as functions 
of 6 (dashed lines). In the latter case, the Routhians are 
drawn by uniformly scaling the range of 0° < < 90° into 
the corresponding ranges of 7. 

Figure [2] shows the results of the SCTAC calculations. 
If the axis of rotation agrees with one of the principal 
axes (PAC), SCTAC coincides essentially with the CNS 
of Ref. (In SCTAC the Strutinsky renormaliza- 

tion is only carried out for the non-rotational part of 
the Routhian, whereas in CNS the rotational energy is 
also renormalized.) The equilibrium deformation param- 
eters in SCTAC for rotation about the short axis are 
£2 = 0.325, 7 = 22° (lower minimum), and for rotation 
about the medium axis £2 = 0.31, 7 = —15° (higher min- 
imum) , which are close to the ones found by means of the 
CNS for the TSDl and TSD2 minima, respectively [Hj]. 

It is clearly seen in Fig. [2] that the PAC minimum at 
7 = —15° becomes a saddle if the rotational axis is al- 
lowed to tilt. The dashed lines show how the energies 
change in a smooth way when tilting the rotational axis 
from short (9 ~ 0°) to medium (9 = 90°) while keep- 
ing £2 and 7 constant. Thus one cannot associate TSD 2 
with the band observed in 158 £r. In addition, the transi- 
tion quadrupole moment for the stable minimum TSDl 
is ^8 eb, which is too small as compared with the exper- 
imental value of ~lleb (see also the discussion in [la]). 

In the case of SCTAC, the tilt angle of the rotational 
axis is defined relative to one of the principal axes of the 
deformed potential in a straightforward way. In SHF- 
TAC calculations, it must be introduced by means of a 
constraint on the orientation of J along with the con- 
straints Im(Q22) = Q2±i = on the orientation of the 
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3: (Color online) Total Routhian surfaces in the 
322) plane for the A-configuration of Table [I] at uj = 



FIG. 

(Q20 

0.7 MeV and 9 = (PAC) obtained with SkM* (left) and 
SLy4L (right) functionals. The corresponding values of 7 are 
marked. 



principal axis of the nucleus defined in terms of the total 
(mass) quadrupolc moment Qiu,- The conditions on the 
corresponding Lagrange multipliers have been derived by 
Kerman and Onishi (KO) [21|. We use the 2D counter- 
part of relation (3.6) of Ref. [2l|, which states that uj 
and J are not parallel (a ^ 0°) if the Routhian is not at 
a stationary point (e.g., PAC). By using the Augmented 
Lagrangian Method of the HFODD code [23|, we have 
checked that the KO conditions are obeyed to a high 
precision for all angles 9. The resulting angles a do not 
exceed 0.1-0.2°, depending on configuration. 



TABLE I: The SHF configurations in 158 Er studied in this 
Letter. Each configuration is described by the number of 
states occupied in the four parity-signature (ir, r) blocks, in 
the convention defined in Ref. [2g| , and also by the transition 
quadrupole moment Q t . The Qt values are from the SHF- 
SkM* calculations at cj = 0.6 MeV and = 0°. 





configuration 


7r r 


Qt (eb) 


A: 


1/(23, 23, 22, 22] ® tt[17, 18, 16, 17] 


- -1 


7.68 


B: 


i/[23, 23, 22, 22] <g> tt[17, 17, 17, 17] 


+ +1 


7.71 


C: 


i/[23, 24, 21, 22] ® tt[17, 18, 16, 17] 


+ +1 


7.36 


D: 


v[23, 23, 22, 22] <g> tt[17, 17, 17, 17] 


+ +1 


10.72 



In this Letter, we studied four different SHF configu- 
rations listed in Table |TJ They all are expected to appear 
near yrast around <320=19 eb. First, in Fig.[3]we show the 
SHFPAC Routhian surfaces for the A-configuration at 
u = 0.7 MeV and 9 = 0°. The Routhians calculated for 
the two functionals, SkM* and SLy4^, are rather similar. 
Clearly visible two minima at 7 = — arctanQ22/Q20 ~ 
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FIG. 4: (Color online) Total Routhians of 158 Er as a func- 
tion of at w=0.5, 0.6, and 0.7 MeV for SkM* (left) and 
SLy4_L (right) functionals. Dashed, solid, dotted, and dot- 
dashed lines represent TSD configurations A, B, C, and D, 
respectively (see Table [I]). Note that the labeling of the con- 
figurations is only valid when 6 = 0° or 90°. Otherwise, 
the signature is no longer a good quantum number. For the 
SLy4i, functional, converged solutions for configurations C 
and D could not be obtained. 



14° and —10° are well separated by a potential barrier of 
- 1 MeV at 7 w 0°. The SHFPAC Routhians are similar 
to those of the SCTAC and CNS calculations of Ref. QJ], 
which have a somewhat smaller barrier of ~0.7MeV. 

In Fig. [4] we allow the rotational axis to tilt (9 ^ 0) by 
starting from SHFPAC solutions with Q22 ~ —4 b. [A ro- 
tation of this shape around the y-axis (9 = 90°) is equiv- 
alent to that of Q22 ~ 4b around the a;-axis (9 = 0°).] It 
can be seen that for the configuration A, the minimum 
that appears in Fig. [3] at Q22 ~ 4b (7 « —10°) is un- 
stable with respect to a reorientation of the rotational 
axis, that is, it represents the saddle-point. On the other 
hand, the lower minimum, at 7 w 14°, remains stable. 
This is consistent with SCTAC calculations of Fig. O A 
similar situation is predicted for configurations C and D. 
At cj=0.5McV, the configuration B has a minimum at 
9 = 90° (7 < 0), but it becomes ^-unstable at higher 
rotational frequencies and a minimum at 9 = 0° devel- 
ops at w=0.7MeV. This interesting behavior, together 
with a discussion of wobbling modes in bands A-D will 
be discussed in detail in a forthcoming paper. 

A detailed search at larger deformations has resulted 
in configuration D having Q20 ~ 29 b, Q22 ~ — 6 b 
(7 = 12°). This structure has similar occupations as 
configuration B, see Table HI but much larger deforma- 
tion. The calculated Qt ~ 10.5 eb well reproduces the 



4 



1 ' 1 ' 


, | , | , | i | , | i | , 
- SLy4 L D - 


- SkM* D 


'- B ...^* f ' 1 
A 




B 







40 50 60 40 50 60 

Angular momentum 



FIG. 5: (Color online) Relative energies of configurations A 
(dashed), B (solid), and C (dotted) with respect to D (dash- 
dotted line at MeV). 



TABLE II: Charge quadrupole moments, transition (charge) 
quadrupole moments, and angular momenta for bands A and 
D calculated with SkM*. 



Band 


u (MeV) 


Qfo (eb) 


-Q22 (eb) 


Qt (eb) 


J(h) 




0.50 


9.01 


1.92 


7.90 


45.13 


A 


0.60 


8.72 


1.99 


7.57 


51.61 




0.70 


8.42 


2.04 


7.24 


57.33 




0.80 


8.09 


2.10 


6.88 


62.84 




0.40 


11.81 


2.30 


10.48 


31.51 


D 


0.50 


12.22 


2.40 


10.83 


48.36 


0.60 


12.13 


2.45 


10.72 


54.82 




0.70 


12.05 


2.48 


10.62 


61.75 




0.80 


11.96 


2.53 


10.50 


71.90 



experimental value of ~11 eb [l6|. We note that band 
D has a larger qua drup ole moment and smaller 7 than 
band TSD3 of Ref. \jjj§. 

Figure [5] compares the energies of configurations A, B, 
and C with that of D. For the SkM* functional (left), 
which has been partly optimized at large deformations, 
it can be seen that at low spins band D lies about 1 MeV 
above band B and 1.5-1.8 MeV above bands A and C. Due 
to its large moment of inertia, band D crosses band B at 
about J ss 55ft and bands A and C at J w 70ft. This 
is not inconsistent with the spin estimates of Ref. (lOj 
based on the feeding analysis of yrast states. While we 
do not expect the relative energies calculated in SHF us- 
ing current functionals to be precise, as evidenced by 
appreciable differences between SkM* and SLy4L predic- 
tions in Fig. O band D is clearly the best candidate 
for the structure observed in experiment. As seen in 



Fig. 21 D-configuration is perfectly stable against the ro- 
tational axis tilting. In Table IH we list the quadrupole 
moments of the lowest stable TSD configuration (A- 
configuration) at positive-7 value and different rotational 
frequencies. Qt is calculated from charge quadrupole mo- 
ments through the relation Q t = + ^J\Q^ Q- 

In summary, we have performed, for the first time, 
TAC calculations within the self-consistent Skyrme 
Hartree Fock model in which the KO conditions for tri- 
axial rotation are strictly obeyed. To address the recent 
puzzling experimental data, we studied the nucleus 158 Er 
at ultrahigh spins. Restricting the direction of the rota- 
tional axis to one of the principal axes of the density dis- 
tribution yields two TSD minima with similar 62 values 
but with positive and negative 7 deformations, similar 
to our SCTAC predictions and the previous CNS calcu- 
lations [IH . Allowing the rotational axis to tilt away from 
the principal axes shows, however, that the higher-energy 
minimum is actually a saddle point; hence, it cannot be 
associated with a physical state. It is the lower-energy 
minimum that represents a TSD band. We have thus 
clarified a long-standing question pertaining to the na- 
ture of positive- and negative-7 bands associated with 
the same intrinsic shape in the PAC approach: the rota- 
tion of a well-deformed, slightly triaxial configuration can 
be either about a short or medium axis, but not about 
both. 

Several TSD configurations differing by proton and 
neutron occupations and quadrupole moments have been 
investigated. In the angular momentum range of 50-70 
ft, they are predicted to have transition quadrupole mo- 
ments of 7-8 eb, which are below the measured values of 
Qt ~ lleb [l6[. We have identified an excited TSD con- 
figuration with a stable positive-7 minimum, which has a 
large transition quadrupole moment of Qt « 10.5 eb that 
agrees well with the experimental value. At spins higher 
than ~ 55ft, this band - different from structures TSD3 
and SD of Ref. [l(| - lies close to the less deformed TSD 
bands, and it is expected to become yrast above J > 70ft. 
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